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ABSTRACT 

Observationally inferred superburst ignition depths are shallower than models predict. We address this dis- 
crepancy by reexamining the superburst trigger mechanism. We first explore the hypothesis of Kuulkers et al. 
that exothermic electron captures trigger superbursts. We find that all electron capture reactions are thermally 
stable in accreting neutron star oceans and thus are not a viable trigger mechanism. Fusion reactions other than 
12 C + 12 C are infeasible as well since the possible reactants either deplete at much shallower depths or have 
prohibitively large Coulomb barriers. Thus we confirm the proposal of Cumming & Bildsten and Strohmayer 
& Brown that 12 C + I2 C triggers superbursts. We then examine the 12 C + 12 C fusion rate. The reaction cross- 
section is experimentally unknown at astrophy sic ally relevant energies, but resonances exist in the 12 C + 12 C 
system throughout the entire measured energy range. Thus it is likely, and in fact has been predicted, that a 
resonance exists near the Gamow peak energy E pk » 1.5 MeV. For such a hypothetical 1 .5 MeV resonance, we 
derive both a fiducial value and upper limit to the resonance strength (u>j)r and find that such a resonance could 
decrease the theoretically predicted superburst ignition depth by up to a factor of 4; in this case, observationally 
inferred superburst ignition depths would accord with model predictions for a range of plausible neutron star 
parameters. Said differently, such a resonance would decrease the temperature required for unstable I2 C igni- 
tion at a column depth 10 12 gem" 2 from 6 x 10 8 K to 5 x 10 8 K. A resonance at 1.5 MeV would not strongly 
affect the ignition density of Type la supernovae, but it would lower the temperature at which 12 C ignites in 
massive post-main-sequence stars. Determining the existence of a strong resonance in the Gamow window 
requires measurements of the 12 C + 12 C cross-section down to a center-of-mass energy near 1 .5 MeV, which is 
within reach of the proposed DUSEL facility. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: neutron — X-rays: bursts 



1. INTRODUCTION 

Superbursts are long, energetic, and rare thermonuclear 
flashes on accreting neutron stars in low-mass X-ray binaries. 
Their durations (~ hours), fluences (~ 10 42 ergs), and recur- 
rence times (~ years) distinguish superbursts from their typi- 
cal hydrogen- and h elium-trigg ered c ounterparts (for reviews, 
see iKuulkersI 12004 ICumming||2005l: IStrohmaver & Bildstenl 
120061) . As of this writin g, astronomers have detected 15 super- 
bursts from 10 sources iKuulkersl^OO^Iin't Zand et al.ll2004t 
iRemillard et"aTJl2005l: lKuulkersll2005riKeek et al.ll2008l and 
references therein). 

The proposal (ICum ming & Bild stenl 120011: 

IStrohmaver &Brownl |2002[) tha t thermally unstab l e 12 C 
fusion dWooslev & Taaml 1 19761: iTaam & Pickruml 119781: 
iBrown & Bildstenl 119 98) triggers superbursts offers a rea- 
sonable explanation of their origin. Cooling model fits 
to superburst l ight curves ( ICumming & Macbeth! l2004t 
Cumming et al.l 120061) as well as obser ved fluences and 
recurrence times (e.g. iKeek et all 120061) suggest ignition 
column depths £j gn « 10 12 gem" 2 , where E = J pdz is the 
radially integrated density. Previo u s superburst ignition 
mode ls (Cum ming & Bildstenl 120011: IStrohmaver & Brown 
2002 1 : ICummingl l2003b iBrownl 12 004: Co oper & Naravan 



2005 



Cumming et al. 2006; Gupta et al. 2007) demonstrated 



that 12 C ignites at £ 
the ocean temperature T 



10 12 gem 2 only if 12 C is abundant and 
6x 10 8 K at that column depth; 
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within existing models of nuclear heating in the neutron 
star crust, such a large temperature requires an inefficient 
neutrino emission mechanism in the neutron star core and a 
low thermal conductivity in the neutron star crust, so that the 
crust is much hotter than the core. 

Recent observations, simulations, and experiments have 
exposed three fundamental problems with this scenario. 
First and foremost is the inference that the ocean is in 
fact too cold for 12 C ignition at the inferred column depth 
Si gn ~ 10 12 gem" 2 . This com es from fits (IShternin et al.1 
120071: |Brown & Cumming 2009) to the quies cent cooling of 
the q u asi-persistent transien t KS 1731-260 dWiinands et al.l 
120021: iRutledge et ail 120021: ICackett et all I2006J), a s ystem 
that also exhibited a superburst dKuulkers et al.ll2002bl) . The 
timescale for the quiescent luminosity to decrease suggests 
that the crust's thermal conductivity is high; as a result, the 
inner crust temperature remains close to that of the core 
even during the accretion outburst. Cac kett et al.l (120081) 
reach the same conclusi on for MXB 1659-29 (see also 
IBrown & Cumm ing 2009 ). In fact, molecular dynamics sim- 
ulatio n results (IHorowitz et al.l2007l 120091: IHorowitz & Berryl 
2009) suggest that the neutron star crust is arranged in a 
regular lattice and therefore has a high the rmal conductiv- 
ity. Neither shear-indu ced viscous heating dPiro & Bildstenl 
20071 IKeek etan i2009) nor deep crustal heating due to elec- 
tron captur es, neutron emissions, a n d pycnonuclear reac- 
tions (e.g. lHaensel & Zdunikl 120081: IHorowitz et"aT1 120081: 
>ta et al.ll2008h can account for the heat necessary to raise 



the ocean temperature t o the required leve l (although see 
iPage & Cummingll2005t [Blaschke et al.ll2008l who consider 
heating in strange stars and hybrid stars, respectively). 

Second, evidence of heavy-ion f usion hindrance at ex- 
treme sub-Coulomb-barrier energies (iJiang etalJl200ll2007l 
2008) implies that the cross-section and thereby the 12 C + 
12 C reaction rate may be orders of magnitude smaller than 
that assumed in the aforementioned superburst ignition mod- 
els. When included in superburst ignition models, heavy- 
ion fusion hindrance increases Ej gn by at least a factor of 2 
dGasques et al.ll2007l) . 

Third, the means by which nuclear burning on the stellar 
surface produces sufficient quantities of 12 C to trigger super- 
bursts is poorly understood. Super burst models require large 
12 C mass f raction s for ignition dCumming & Bildstenl 120011: 
Cummind 120031: IC ooper & Naravanl 120051: ICooper et al. 
2006t ICumming et al.1 120061) . All systems that exhibit su- 



perbur sts show helium-tri ggered type I X-ray bursts as well 
(e.g. |Galloway et al. 2008), but theoretical models of such 
bursts yield 12 C ma ss fractions far smaller than t h ose re- 
quired for ignition dlossl 119781: ISchatz et al.l [20011 l"2 003b; 
Koike et al.1 |2004| IWooslev et al.l l2004at iFisker et al.l l2005l 
20081: iPeng et all 12007b iParikh et al.l 120081) . Most sys- 



tems that exhibit superbursts apparently undergo long pe 
riods of stable nuclear burning between successive helium 
triggered bursts (Kuulkers et al. 20 02atlin 't Zand et alJl2003t 
Keek et al. 2008); stable burning generates much more 12 C 
than unstable burning, but the calculated yield is insufficient 
to trigger superbu rsts in all systems, particularly those accret- 
ing at a high rate dTaam &Pickluril 1978b ISchatz et al.lll999l 
I2003bt ICooper et al.ll2006HFisker et al.ll2006l) . 

Detection of a superburst from t he classical tra nsient 
4U 1608-522 jRem illard et all [2001 iKuulkersI 120051) with 
Sign ~ 10 12 gem" 2 (Keek et all 120081) exacerbates all three 
problems: (1) the transient's inferred ocean temperature is 
lower than those of other systems exhibiting superbursts, (2) 
heavy ion fusion hindrance is greater at lower temperatures, 
and (3) most of the matter accreted onto the neutron star 
prior to the observed superburst likely burned during helium- 
triggered type I X-ray bursts, which current theoretical calcu- 
lations suggest would generate far less 12 C than that required 
for ignition. 

Reconciling superburst observations with the current the- 
oretical model is impossible. This motivates both a critical 
assessment of the current ignition model and a search for al- 
te rnative ignition mechan isms. 

iKuulkers et alJ (2002b) proposed such an alternative mech- 
anism. Electron captures onto protons and the subsequent 
captures of th e resulting neutrons onto he avy nuclei liberate 
« 7 MeV/m u dBildsten & Cummingll998l) . Prethreshold cap- 
tures of super-Fermi electrons are very temperature sensitive 
and therefore could trigger an energetic thermonuclear flash. 
An attractive feature of this mechanism is that ignition oc- 
curs always at the same electron chemical potential; thus E; gn 
would be similar for all superbursts, in accord with observa- 
tions. Unfortunately, the calculated Ej gn «2x 10 10 gem" 2 
is much smaller than the inferred superburst E; gn , making 
it an unlikely trigger mechanism. This motivates our inves- 
tigation of exothermic electron captures onto heavy nuclei, 
which occur throughout the ocean and crust of an accr eting 
neutron star, including the superburst ignition region (ISatol 
1979tlBlaes et aLlll990tlHaensel & ZdunikTl9 90. 2003,|200! 
Gupta et al.l2007l) . In ^2] we determine the thermal stability of 
electron captures onto heavy nuclei in accreting neutron stars. 



We show that instability requires unrealistically large reaction 
Q-values, where Q is the energy released per capture; thus 
we conclude that electron captures in accreting neutron star 
oceans are thermally stable. We then consider the relevance of 
a captures onto light elements such as 12 C in §|3] We find that 
none of these reactions is a feasible mechanism and thereby 
confirm the proposal that 12 C + I2 C triggers superbursts. 

In 50]we assess whether the 12 C+ 12 C reaction rate could be 
much larger than the fiducial rate. We investigate (§ 14. It the 
screening enhancement factor, including a careful evaluation 
of corrections to the liner mixing rule, and show that uncer- 
tainties in the plasma screening enhancement are unlikely to 
change the 12 C+ 12 C reaction rate enough. We then consider 
the nuclear cross-section. We find that a strong resonance at 
an energy near 1.5 MeV in the I2 C + 12 C system, which theo- 
retical nuclear physics models predict, could increase the re- 
action rate in the astrophysically relevant temperature range 
by over two orders of magnitude. In §|5]we show that the ex- 
istence of such a resonance could decrease the predicted Si gn 
by a factor « 2—4 and thereby alleviate the discrepancy be- 
tween superburst models and observations. We conclude in 
$6]by discussing the implications of our findings. 

2. THERMAL STABILITY OF ELECTRON CAPTURES 

Consider the accretion-driven compression of a matter ele- 
ment containing a nucleus of mass M(A,Z), where A is the 
mass number and Z is the proton number. The degener- 
ate electrons' chemical potential rises as the nucleus ad- 
vects to higher pressures. Eventually M(A,Z)+fi e /c 2 exceeds 
M(A,Z- 1) and electron capture becomes energetically fa- 
vorable. Such captures often occur in equilibrium and re- 
lease a negligible amount of energy; however, some cap- 
tures can heat the ocean in two mutually inclusive ways (e.g. 
Gupta et al. 2007): (1) An electron captures into an exited 
state of the daughter nucleus if, for example, the daughter nu- 
cleus's ground state is forbidden. The daughter nucleus then 
radiatively deexcites and thereby heats the ocean. (2) If the 
parent nucleus is even-even, then M(A,Z- 1) > M(A,Z-2) 
due to the nuclear pairing energy, and a second electron cap- 
ture immediately ensues. The latter, post-threshold electron 
capture occurs out of equilibrium and thus releases heat. 

2.1. Governing Equations 

We construct a simple model of the accreted layer to deter- 
mine the stability of exothermic electron captures to thermal 
perturbations. We assume spherical accretion onto a neutron 
star of mass M = 1.4 M Q and radius R = 10 km at an accre- 
tion rate per unit area S. The accreted layer's scale height 
is much less than R, so we set the gravitational acceleration 
g = GM/R 2 (l-2GM/Rc 2 T^ 2 = 2.43 x 10 14 cms" 2 through- 
out the layer. The layer is always in hydrostatic equilibrium, 
so the column depth E is a good Eulerian coordinate. To fa- 
cilitate comparisons between microphysical and observation- 
ally inferred quantities, we express microphysical quantities 
in terms of the macroscopic coordinate £ using the follow- 
ing approximate relation between mass density p and E for 
relativistic, degenerate electrons 1 , 



P' 



:5.9xl0Vm-(^W 4 , 



(1) 



1 In this and following expressions, we suppress the scaling with g and 
evaluate the expressions at g = 2.43 X 10 14 cms -2 . 
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where (A/Z) is the mean molecular weight per electron and 
E = E12 x 10 12 gem" 2 . We denote the Eulerian time and 
spatial derivatives as d/dt and d/dY,, respectively, and the 
Lagrangian derivative following a matter element as D/Dt, 
where D/Dt = d/dt + E<9/<9E. The governing transport, en- 
tropy, and continuity equations are 



dT 

F = pK ^ 

Ds „ 
T- = 6Xr a 

DX 

~b7 



dF 

5E : 



= -Xr, t 



(2) 



(3) 
(4) 



where F is the flux, K is the thermal conductivity, s is the en- 
tropy, £ = Q/(Am u ) is the energy per gram released via elec- 
tron captures, X is the parent nucleus mass fraction, 



In 2 

W) 



1 



(m e c 2 ) 5 J Q 



e\e- 



Qff(E,^,T)dE (5) 



is the electron capture rate (Full er et all 1 19851 ), (// 
is the effective ft value dFuller et all 1 19801 1 19851 
Langa nke & Martmez-Pinedol 1200 ll) . m e is the electron 
mass, Q is the threshold energy, and 



f(E,fi e ,T) = 



1 



(6) 



l+exp[(E-fi e )/k B T] 

is the Fermi-Dirac distribution function. 

Consider prethreshold electron captures, where p, e < Q. For 
T = 0, all electrons have energies E < fi e by equation (O; 
electron capture is blocked. For T > some electrons have 
E > Q and thus can capture. The number of electrons with 
E > Q increases with T, which makes prethreshold electron 
capture temperature-sensitive. In the prethreshold limit (p e - 
Q)/k B T < 0, 

2Q 2 (k B Tf 



, ln2 

r,= l w 



-exp 



(7) 



(m e c 2 ) 5 

(iFuller et al.lll985t iBildsten & Cumminglll998l) . Conversely, 
for fi e > Q a majority of electrons has E > Q and hence 
can capture for any T, making r ec relatively temperature- 
insensitive and thus thermally stable. We therefore consider 
prethreshold electron captures exclusively hereafter. 

2.2. Prethreshold Electron Capture 

Prethreshold electron captures occur within a thin layer 
in the deep ocean. To illustrate this, consider the height- 
integrated capture rate. Relativistic, degenerate electrons sup- 
ply the pressure P = gE, so 

,4 



= 10 12 gcm- 2 ( ; 



Me 



(8) 



12n 2 g(hc) 3 " °'~ V3.43MeV< 
Integrating equation (0 over E and using equation (H), 

r ec (^,T)dE' w r ec (/i e ,T)E (9) 

V Me / 

for /i e 3> k B T. Equation (O shows that prethreshold electron 
captures occur in a narrow column depth range 



AE 4k B T 



E /Lie 

where T = n x 10 8 K. 



= 0.050 



Tg\ 1/4 



J 12 ! 



(10) 



Now consider electron captures in steady-state, such 
that electrons capture onto nuclei at the same rate as 
accretion advects the nuc lei (see also the discussion in 
IBildsten & Cummindll998l) . Equation (0]i becomes 



■ dlnX 
L ^E~ 



(11) 



Integrating equation ( fTTT ) from to Q and using equations I© 
and ©, we find that most electron captures occur prethreshold 
when 



E12 > 0.027 



-16/5 



4/5 



0.3E E dd J 



(ft) 
10 3 s 



4/5 



(12) 



where Eedd ~ 10 5 gcm" 2 s _1 is the local accretion rate at 
which the accretion flux equals the Eddington flux 2 . Super- 
bursts ignite at column depths E12 ~ 1 and accretion rates 
E w 0.1-1 EEdd- Equation ( fT2l shows that superallowed elec- 
tron captures (for which (ft) ~ 10 3 — 10 4 s) occur prethreshold 
at superburst ignition depths. 

2.3. Thermal Stability Analysis 

We now derive a o n e-zone model (e .g., Fuiimot o" et alj 
11981b iPaczynskil Il983t IBildsten! Il998bl) from the govern- 
ing equations (HJU) to determine the stability of prethreshold 
electron captures to thermal perturbations and thereby ascer- 
tain whether electron captures trigger superbursts. We con- 
sider only temperature perturbations and ignore the accretion- 
induced entropy advection through the bottom of the zone. 
Therefore, we set d/dt = in equation (01 and approximate 
Ds/Dt = ds/dt in equation ©. Perturbations occur at con- 
stant pressure since the scale height ~ E/p <C R; therefore, 
we write Tds = CpdT, where Cp is the specific heat at con- 
stant pressure. Equations (f2]|4| become 



F = pK 



,8T 

as' 



dT dF 
C P —-£Xr ec + — : 

■ dX_ 

9E 



(13) 



(14) 



(15) 



We simplify equations ( TT~3TT~5b as follows. We set p, F, and 
X to be constant throughout the layer; specifically, we adopt 
step-like profiles for F and X: 

F(E,f) = J F (f)e(E ec -E), X(E)=Z e(E ec -E), (16) 

where O is the Heaviside step function, E ec denotes the col- 
umn depth at the bottom of the layer, and Fo and Xq denote 
the values at the top of the layer, where E <C E ec - We as- 
sume the ocean consists of a single ion species and set Xq = 1 . 
Ele ctron-ion scattering sets the ocean's thermal conductiv- 
ity jYakovlev & Urpin|[T980t lltoh et al] [19831: IPotekhin et alj 
ri999h 



9.8 x 10 18 



p 



6 x 10 s gem" 3 



1/3 



ergs cm 1 s K , 
(17) 



2 We define the Eddington luminosity to be 47rGMc/«es in the frame of a 
dist ant observer. Thi s is the largest luminosity observable by such an observer 
(see Shapiro & Teukolsky 1983). This differs from the definition used by 
IGallowav et aTl J2008|) by a factor of [1 -2GM/(Rc 2 )T 1 / 2 = 1.3. 



4 



where we set the Coulomb logarithm A e ; = 1, a value appro- 
priate for a plasma at T% rj 5 and p w 6 x 10 s gem -3 . Since 
oc 7\ we rewrite equation ([T3) as 

_ pK dT 2 



2T <9E ' 

Equations (16) and ( fl~8T > then imply 

r 2 (s,o = 7; ) 2 +[7; c (f) 2 -r () 2 ] 



(18) 



(19) 



Integrating equations (1141115b over E and using equations ©, 
CO), (H6), (ED, and (Hi, we find 



C, 



57; c AE 



ft 



E e . 

£ 

AE' 



£r e . 



pK 
2T 



-T 2 



(20) 



(21) 



Equation (2T) shows that electron capture occurs when the 
lifetime of a parent nucleus l/r ec equals the time f acc AE/E 
an advecting element spends within the capture region, where 

r acc = £/£ (22) 

is the accretion timesc ale. Note that t his diffe r s from 
the usual assumption (iBlaes et all Il990t iBildstenl Il998at 
iBildsten & Cummingl Il998t lUshomirskv et all l2000h that 
electron capture occurs when l/r ec = f acc . 

Finally, conducting a linear stability analysis on equation 
(20) and using equation (21), the thermal instability criterion 
is 



v — > 

fa, 



S5L 



(23) 



where the temperature-sensitivity of the height-integrated 
electron capture rate 



dlnOecAE) 
d\nT 



:4 + ln 



81n2 



Q 2 (k B T) 4 



(ft) (m e c 2 )Ve 



(24) 



from equations (7) and (9). Noting that electrons capture 
when p s /Q~ 1, we write equation (24) as 



:8.14 + ln 




( ift) 



V 10 3 : 



(25) 

Using equations (22) and (fT7) , the thermal instability criterion 
(23) becomes 



>20 



.V 



0.3EEdd / 



A 
2Z 



MeV, (26) 



where Q is the energy released per electron capture. 

We tested the accuracy of equation (23 > using the suitably 
modifi ed global linear stability analysis of lCooper & Naravanl 
(20051) . The minimum Q for instability derived from the 
global stability analysis differed from that of equation (26) 
by less than 30% for each of the 12 test cases. 

Typically Q < Q because the daughter nucleus is generally 
more massive than the p arent nucleus, although exceptions 
exist. iGupta et all (120071) find Q < 6.2 MeV for all electron 
captures that occur for p e < 6 MeV, or equivalently, Era < 10 
(eq. (8)). From equation (26) . it follows that electron captures 
are thermally stable for the accretion rates and column depths 
at which superbursts occur. Therefore, we conclude that elec- 
tron captures do not trigger superbursts. 



3. ALTERNATIVE FUSION REACTIONS 

In this section, we examine whether light-element fusion 
reactions trigger superbursts. Hydrogen has an electron cap- 
ture threshold energy Q = 1.2933 MeV and thus depletes at 
E12 < 2 x 10" 2 (eq. 0). The helium abundance at £jg n is 
less certain, and we discuss it below. The paucity of stable 
isotopes of Z = 3-5 nuclei leaves I2 C as the next reasonable 
alternative, which we address in §|4] Finally, nuclei with Z > 6 
are unlikely candidates because the extra Coulomb repulsion 
causes the fusion rates to be significantly lower than that of 
12 C. 

Thus, other than 12 C + 12 C, a capture reactions such as 
12 C(a, 7) ie O are the only plausible fusion reactions that might 
trigger superbursts. The conditions for these reactions to pro- 
duce superbursts are similar to those for electron captures, 
namely, that (1) the reaction rate r nuc is sufficiently temper- 
ature dependent to produce unstable burning, and (2) a par- 
ticles must survive to the inferred Ei gn . Below, we show that 
the latter condition is not met; thus a particles deplete too 
quickly to trigger superbursts. 

The condition for a particles to survive at a column depth 
E is Y/J2i( r mc)i > E/E (see gO) , where Y is the helium 
mass fraction and the sum is over all reaction s that consume a 
particl es. Using the triple-a reaction rate of Fushiki & Lambl 
(19871) and setting p = 6 x 10 8 gem" 3 , T s = 5, and Y = 1, we 
find r nuc = 2.2 x 10 6 s , which implies a lifetime of 4.6 x 
10~ 7 s. The reaction rate r nuc oc T 3 , so the lifetime is much 
larger for smaller helium abundances. The accretion timescale 
f acc = 10 7 [E 12 /(E/E Edd )] s (eq. E2), indicating that Y < 10" 7 
for helium to survive. At this low Y, the rise in temperature 
from consuming the helium via, e.g., 12 C(a, 7), is < 10 6 K <C 
T and hence insufficient to trigger a thermal instability. 

From the results of this section and $2] we conclude that 
12 C fusion triggers superbursts. 

4. THE 12 C + 12 C REACTION RATE 

A possible solution to the superburst ignition problem is 
that the true 12 C + 12 C fusion rate is larger than assumed. 
12 C ignition at the inferred E; gn requires a ~ 10 4 reaction 
rate enhancement for an ocean temperature of 4 x 10 s K 
(Cumming et al.l2006l) or a ~ 10 2 enhancement for 5 x 10 8 K. 
The two sources of uncertainty in the fusion rate are (1) 
plasma screening effects and (2) the nuclear cross-section 
a(E). In the following subsections, we investigate whether 
either source could account for such a large increase in the 
fusion rate. 



4.1. Plasma Screening 

Superbursts ignite in a strongly-coupled Coulomb plasma. 
Two dimensionless parameters determine the plasma's state. 
The first is the Coulomb coupling parameter 



r : 



aksT V 5 



.1/4 



5/3 



(27) 



where a = (3Z/47rn e ,) 1//3 is the ion-sphere radius, n e is the elec- 
tron number density, and we used equation (T), which as- 
sumes the gravitational acceleration g = 2.43 x 10 14 cms -2 . 
For r< 1 Coulomb coupling is weak and the ions constitute 
a Maxwell-Boltzmann gas. As T increases, the ions gradually 
beco me a Coulomb liquid. Wh en T > 175 the ions crystal- 
lize (Potekhin & ChabrierjfeOOOl) . Equation (27) implies that 



5 



superbursts ignite in a Coulomb liquid. The second dimen- 
sionless parameter, 



3r 



; 0.17 



-2/3 



4/4 
J 12 



2Z\ 



1/3 



(28) 



r \ 5 

is the ratio of the classical turning point to the ion separation, 
where 

r=| w f. ^ " (29) 



2£ B 7^ / V 4A; B 7 , /j / 
and /i is the reduced mass of the reacting nuclei. Specifi- 
cally, £ = rjp /a, where rxp is the radius at which the Coulomb 
energy Z 2 e 2 /ry p equals the classical Gamow peak energy 
dClavtonll 19831) 

rk B T 

(30) 



£P k = 



Many-body interactions in a strongly-coupled Coulomb 
plasma modify the Coul omb potential b etween two reacting 
nuclei (for a review, see Ic himar u 1993). From these many- 
body interactions, one derives an effective two-body potential 

zV 

V(r)= H(r), (31) 



where r is the distance between the reacting nuclei. From 
H(r), one derives the plasma screening enhancement to the re- 
action rate exp{(H{r)) /k B T), where ( H(r)) is a path-integral 
average of H(r) (e.g. Ilchimarul [19931) . One can expand the 
static mean-fie ld potential H{r) as a power series in (r/a) 2 
dWidomll 1963b . Neglecting quantum effects in H(r), the lead- 
ing order term H(Q) is a thermodynamic quantity; H(0) equals 
the difference between the Coulomb (o r excess) Helmholt z 
free energy before and after the reaction dDeWitt et al.ll 1 973b . 

Monte Carlo simulations and hy pernetted chain calcu- 
lation s of binary ionic m ix tures dHansen & V ieillefosse 
19761 lHansen et all fl977l IChabrier & Ashcrofj 11990 



Ogata et al.ll 1 9931: [Rosenfeldll 1 9951 1 1 996t iDeWitt et alii 1 996 ; 
DeWitt & Slattery|2003l) suggest the excess free energy obeys 



the linear mixing rule to high accuracy in the regime F > 1 
dPotekhin et al.1 120091) . Therefore, authors usually invoke 
the linear mixing rule when deriving the plasma screening 
enhancement to the reaction rate. In this case, the total free 
energy of an ionic mixture 

z^ex 

r ^ = V^«r x (r,), (32) 

k B T ^ 

I 

where N, is the number of ions with charge Z,, / ex = 
F ex ' 0CF /Nk B T is the well-determined reduced excess 
free energy per ion of a one-component plasma (e.g., 
Chabrier & Pot ekhinl [T998t iPotekhin & Chabrieril2000l) . and 

T[ oc Zj is the Coulomb coupling parameter for species i 
(eq. E3). 

From equatio n ([32j, H(0)/k B T = 2f ex (T) - f ex (2 5 ^T) 
(lJancovicill977l see also the Appendix), where T is that of the 
reacting ions; using the ion-sphere model result / ex w -0.9r 
dSalpeterl 19541) . H(0)/k B T = 0.9(2 5 / 3 -2)r « 1.0573r, so the 
lowest-order screening enhancement to the I2 C + 12 C reaction 
rate 

-l 



exp 



(H(0)\ = 



5.5 x 10 2 exp 



4/4 



(33) 



Despite its simpli city, equation <f33T> is adequate for 
most applications dGasques et al. 2005; Yak ovlev et al.ll2O06i; 
Chug unov et al.ll2007l) . We discuss corrections to the screen- 
ing enhancement below. 



4.1.1. Deviation from Linear Mixing Rule 

The excess free energy F ex of a multicomponent plasma ex- 
hibits small deviations from the linear mixing rule (eq. | [32l ). 
In general, 



ks~f 



Nif^Td+NAf 



(34) 



where A/ ex > is a function of both the charges Z, 
and concentration s jc, = N,/N of the ionic species (e.g., 
IDeWitt eTa l. 1996). Using the hypernetted chain calculations 
oflDeWitt et all dl9£6j) an d the ansatz A/ ex oc x l x 2 (Z 2 /Z l ) i l 2 
(DeWitt & Slattery 2003), we find 



A/ ex = (0.0091 lnrj+O.OlS^^ I ^ 

. Zi 



3/2 



(35) 



for a binary ionic mixture (see also Potekhin et al. 20091). 

To incorporate linear mixing rule deviations into H(Q) cal- 
culations, previous authors assumed a one-component plasma 
(consisting of 12 C ions in this case). Fusion of two 12 C ions 
generates a compound nucleus 24 Mg and thereby forms a bi- 
nary ionic mixture. One determines H(0) by finding the dif- 
ference in F ex before and after the reaction in the limit tha t 
the compound nucleus concentration x 2 — > (Ichimaru 1993). 
Using this assumption and equations d34l ) and d35l ). we find 
the correction to H(0): 

AH(0) 



k n T 



= -0.026 lnr-0.051, 



(36) 



in good agreement with IDeWitt et al.1 d 1996b ; for a one- 
component plasma, linear mixing rule deviations reduce the 
plasma screening enhancement factor by ~ 10%. 

However, the plasma at the superburst ignition depth is 
likely a mixture of 12 C and heavier ions with Z < 46 
(e.g.. lKoikeetal.1 fl999l 12001; ISchatz et all l200ll 120031 
IWooslev et al.l2004al) . Generalizing equation d35T l for a multi- 
comp onent plasma with Z, < Zj for i < j, we find dOgata et al.1 
[19931) 



3/2 



A.r = ^uA/:;: A/^ = (0.00911nr,-+0.018) M 

i<j ^ ' 

(37) 

To illustrate the effect spectator ions have on the screening 
enhancement, consider for simplicity a ternary ionic mixture 
of 12 C, 24 Mg, and a representative spectator ion 56 Fe. Using 
equations d34T > and d37l > and again taking the 24 Mg concentra- 
tion X2 — » 0, so that x\ +X3 = 1, we find (see Appendix) 

=-A/f x +*3 [A/r x + (l+x 3 )A/--A/ 2 ex ] . (38) 

Note that equation d38l reduces to d36] l in the limit Xj — > 0, 
as it should. Equation ( f3~8b shows that, since the bracketed 
term is positive, heavy spectator ions increase the screening 
enhancement factor (i.e. dAH{0)/dxi > 0). For the fiducial 
12 C mass fraction 0.2 a nd 56 Fe mass fraction 0.8 (such that 
X\ =7/13 and JC3 = 6 /13; ICumm ing et al. 2006), linear mixing 
rule deviations increase the plasma screening enhancement by 

» 10%. 

Pote khin et al.1 (120091) developed an analytic formula for 
A/ ex that is more accurate than equation d3~71 >. We derive the 
corresponding formula for AH(Q)/k B T in the Appendix. In 
Figure Q] we plot AH(0)/k B T as a function of the specta- 
tor ion number fraction X3 using both equation d38l and the 
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*3 

FIG. 1. — Correction to the plasma screening enhancement factor due to 
linear mixing rule deviations as a function of the spectator ion number frac- 
tion *3. Considered is a ternary ionic mixture with Z\ = 6, Zi — 12, Z3 = 26, 
and Ti = 6. The number fraction of the product ion X2 = 0, so X\ = \—x?>. 
"PCR0 9" refers to the expression for AH(0)/k^T derived from the results of 
IPotekhin et al. 12009, see Appendix). 



expression derived from lPotekhin et all (120091) . Figure[T]con- 
firms that heavy spectator ions increase the plasma screening 
enhancement to the reaction rate, although the two expres- 
sions for AH(0)/kftT differ quantitatively. 



4.1.2. Corrections to (H(r)) 

The next term in the expansion of H(r) goes as (r/a) 2 oc 
C 2 ; its contributi on is small because ( 2 <C 1 (eq. 
iJancovicild 19771) . we find 



). From 



(H(r))-H(0) 



32 



rc 2 



-0.027 



-7/3 



.3/4 
J 12 ■ 



(39) 



This result agrees very well with more accura t e calculation s 
dAlastuev & Jancovic|[l978l: lOgata et aljfl99ll IOgatalll997l) . 
Higher order terms are even smaller; therefore, we conclude 
that corrections to (H(r)) are unimportant in calculating the 
plasma screening enhancement for 12 C ignition. 

4.1.3. Electron Screening Corrections 

In the above analysis, we tacitly assumed a uniform 
electron density. Although highly degenerate, electrons 
nonetheless slightly concentrate around positively charged 
ions. Electron polarization mitigates the Coulomb repulsion 
between ions relative to an unpolarized configuration. This 
has two counteracting effects: It (1) lowers the Coulomb 
repulsion between the two reacting ions, which increases the 
reaction rate, and (2) attenuates the many-body Coulomb 
interactions and thereby H(r), which decreases the reaction 
rate. The Yukawa potential Z 2 e 2 /rexp(-r/rxF) describes 
the two-body potential, where rjp is the Thomas-Fermi 
screening length. For relativistic, degenerate electrons, 
r TF /a = 3.0(Z/6) -1 / 3 (e.g., lHaensel et aj1l2007l). so electron 
screen ing is weak; from the results of Sahrling & Chabrierl 
( 1998), electron screening changes the reaction rate by < 1%. 

Corrections to the lowest-order plasma screening enhance- 
ment (eq. 11331 ) change the 12 C + C reaction rate by a factor 
< 2. Therefore, we conclude that uncertainties in the plasma 
screening enhancement are too small to explain the discrep- 
ancy between superburst observations and theoretical model 
results. 



4.2. The Nuclear Cross-Section 

Although the plasma screening enhancement to the I2 C + 
12 C reaction rate is well-determined for superburst conditions, 
the nuclear cross-section <j{E) is not. Many groups have 
measured a{E) at various center-of-mass energies E down 
to k, 2.1 MeV dPatterson et all 1 19691: iMazarakis & Stepheni 
1971 [1971 LSninka & Winkled fl97i Iffigh & Cuied \WTi 



Kettne retall fl977l 119801: lErbetal.1 Il980t iTreuetalJ 
1980t iBeckeretalJ 119811; iDa smahapat ra et alj 119821; 
Satkowiak et al. 1982; Rosales et al. 2003; Barron-Palos et al. 
2004 120061; lAguilera et al.ll2006t ISpillane et alj|2007l) . How- 



ever, the energy range of interest is centered at the classical 
Gamow peak energy (cf. eq. l29l - [30l ) 



E pk = 1.5 



2/3 



MeV 



and has a full width 



A£ pk = 4 



/ E pk kftT 



1/2 



= 0.59 



5/6 



MeV. 



(40) 



(41) 



Thus, a(E) in the astrophysically relevant energy range is ex- 
perimentally unknown. 

This situation is common in nuclear astrophysics: to deter- 
mine the astrophysical reaction rate, one either extrapolates 
the experimental dat a to lower energies or calculates the rate 
theor etically (e.g., Caug hlan & Fowler! 119881 ; iGasques et alj 
120051) . In doing so, one tacitly assumes the astrophysical rate 
has the nonresonant form, i.e., no prominent resonances exist 
in the compound nucleus within the relevant energy range 3 . 
However, several groups have detected strong resonances in 
the 12 C + 12 C system at energies below the Coulomb bar- 
rier. Resonances exist throughout the entire energy range 
probed so far, and the spacing between adjacent resonances 4 
is ~ 0.3 MeV. Therefore, a resonance probably exists near 
£ pk dBromlev et al.lll960t l Almqvist et aljfl960l: iGalster et alj 
19771: iKorotkv et al.l ll979i lErb et alj|1980t iTreu et alj|1980h 
Spillane et all 120071). In deed. iMichaud & Vogtl d!972|) and 
Perez-Torres et al.f i[2b06) predict a resonance exists in the 12 C 
+ 12 C system with energy £r « 1.5 MeV. If the resonance is 
strong, the thermally averaged reaction rate (av) would be 
much larger than assumed. 

To illustrate the effect a strong resonance within the Gamow 
window would ha ve on the 12 C + 12 C reactio n rate, we follow 
the prediction of iPerez-Torres et"ail (|2006) and assume the 
existence of a single, narrow resonance with Er= 1.5 MeV. 
Then 



(uv) = (<7V)nr+(<7V)r, 



(ctv)r = 



3m u fc B r 



3/2 



h (w7) R exp 



E R 
~k B T 



(42) 



(43) 



3 This statement is not strictly true for heavy ion fusion reactions such 
as 12 C + 12 C. The compound nucleus 24 Mg has numerous quas i-stationary 
states at excitation energies near 15 MeV above the ground state <Endll990l : 
IFirestone1l2007l) . where E pk lies; thus all reactions are resonant. However, 
when the mean level spacing of quasi-stationary states D <k$T, one com- 
putes an average cross-section o ver all resonances, and the reaction rate as- 
sumes the nonresonant form (e.g. Cameron 1959; Fowler & Hovle 1964) 

4 The average level spacing of the detected resonances is much greater 
than that of the quasi-stationary states in the compound nucleus 24 Mg. Thus 
the observed resonances are not ordinary compound nuclear states. As 
lAlmqvis t et al. 1 1 960) first suggested, the resonances are probably quasi- 
molecular doorway states in the 12 C + 12 C system (e.g., Betts & Wuosmaa 

EH. 
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where (<jv)nr is the nonre sonant contribution to the t otal re- 
action rate as given in, e.g., ICaughlan & Fowled (I 19881) . (av) R 
is the resonant contribution, 



(w 7 ) R = 2(27+1) 



r c (r R -r c ) 



' 2(27+ i)r c 



(44) 



is the resonance strength, 7 is the total angular momentum of 
the resonance, Tc is the entrance channel width, and Tr 3> Tc 
is the resonance width. 

The resonant contribution (ctv)r oc (uj^)r. Using the Breit- 
Wigner single resonance formula, 



a(E) = 



ttH 1 



(w7) R r R 



12£m u (£-£R)2 + (r R /2)2 

(e.g.. lClavtonlll983l) . Evaluating equation (@5j at E = E R , 

« / Tr \ / <r(E R ) \ / E R 
(w7)r = 3.4x10 -8 ' 1 ' w 



(45) 



100 keV 



10" 13 barn 



1.5 MeV 



eV, 



(46) 

where we normalize the resonance width T R to that typical 
of kn own resonances (see, e.g., Table IV of lAguilera et al.l 
2006) and the c ross-section at res o nance a(E R ) to the ap- 
proximate value iPerez-Torres et al.l (120061) predict. For this 
work, we adopt (u)j) R = 3.4 x 10~ 8 eV as the fiducial reso- 
nance strength. 

To determine an upper limit for (cJ7) R , we demand that 
the resonance's contribution to the astrophysical S-factor at 
a given energy E, S R (E), be less than the experimentally mea- 
sured value S exp (E) for all E > 2.1 MeV, the lowest energy 
probed at the time of this writing. The 5-factor for 12 C + 12 C 
is 



S(E) = a(E)Eexp 



87.21 



e y 1/2 / 

j +0.46' 



MeV 



VMev) 



(47) 

(Patte rson et~aT1ll969b IClavtonlll983l) . From equations (05) 
and (4T\ . we write 



S R {E) = S(E R )- 



(r R /2) 2 



'(£-£ R ) 2 +(r R /2)2- (48) 

Using equations (gU, (07)1, and (gHJ, demanding that S R (E) < 
S exp (E), and noting that (E-E R ) 2 > (r R /2) 2 , we find 

r R 



(w 7 ) R <5.5 x 10" 



E-E R 
MeV 



100 keV 

2 



-1 



SexpOE) \ 



10 16 MeVbarn 



eV (49) 



for a resonance at E R = 1 .5 MeV. Equation d49l must be sat- 
isfied for all E. According to the experimental data, the min- 
imum value of the b racketed term is r* 1 (see, e.g., Fig. 4 of 
Snill ane et al]l2007h . so 



(u<y) R < 5.5 x 10" 



r R 



100 keV 



eV. 



(50) 



If r R w 100 keV, then our fiducial strength (ujj) r = 3.4 x 
10~ 8 eV is comparable to the maximum possible strength. T R 
may be much smaller, however; the resonance at 2.14 MeV, 
the lowest-energy re sonance known as o f this writing, has a 
width T R < 12 keV dSpillane et al.ll2007l) . Therefore, we set 
(co>7)r = 3.4 x 10~ 7 eV, which is ten times larger than our fidu- 
cial rate, as a reasonable upper limit. 




FIG. 2. — Ratio of the total thermally averaged reaction rate (<xv) = 
(ot>)nr + {&v)r to the nonresonant contribution (<tv)nr for a hypothetical 
1.5 MeV resonance with strength (li?7)r = 3.4 X 10~ 8 eV, the fiducial value, 
as a function of temperature T = Tg X 10 s K. The resonance increases (trv) 
by a factor > 25 near T% 5. 

Figure |2] shows the effect a 1.5 MeV resonance has on the 
reaction rate (trv). For the fiducial (luj) r value, the resonance 
increases (av) by a factor > 25 at temperatures relevant to su- 
perbursts; for the (o>7) R upper limit, the resonance increases 
(av) by a factor > 250. These increases are of the order re- 
quired to reconcile the observationally inferred S; gn with that 
calculated from theoretical models for a specific range of as- 
sumed crust thermal conductivities and core neutrino emissiv- 
ities. In the following section, we compute the superburst S; gn 
with the effect of this resonance. 

5. EFFECTS OF A RESONANCE ON SUPERBURST IGNITION 

We use the global linear stability analysis of 
Cooper & Naravan (2005) to determine the effect a strong 
resonance in the 12 C + 12 C system would have on the 
superburst ignition depth £j gn . We assume steady spherical 
accretion onto a neutron star of mass M = 1 .4 M and radius 
R = 10 km. The accreted matter composition is that of the 
Sun: The hydrogen mass fraction X = 0.7, helium mass 
fraction Y = 0.28, and he avy-element ma s s fract ion Z = 0.02. 
Furthermore, we follow Cummin g et al.l d2006l) and assume 
the 12 C mass fraction Xq = 0.2 at the base of the accreted 
layer. 

We make the following two modifications to the mod el of 
ICooper & Naravanld2 005). (D lCooper & Naravanld2005l) fol- 
lowed Brown (2000) and assumed the energy generated by 
electron captures, neutron emissions, and pycnonuclear reac- 
tions in the crust was distributed uniformly between £12 = 
6 x 10 3 and 2 x 10 5 . We now follow lHaensel & Zdunikl d2008l) 
and distribute the energy according to their Table A. 3. (2) 
Plasma screening reduces the entrance channel width Tc- 
Therefore, the plasma screening enhancement for the resonant 
contribution to the reaction rate includ es a correction factor 
that reduces the ov erall enhancement (Sal peter & Van Hornl 
119691: iMiflerl 119771) . although the reduction is only a few 
percent for the conditions rel evant for superbursts (see, e.g., 
Fig. 1 of ICussons et alJl2002h . We now use the formalism of 
lltoh et all ([2003) for the plasma enhancement factors of both 
the resonant and nonresonant contributions. 

The 12 C + 12 C reaction rate, accretion rate E, and ocean 
temperature profile together determine E; gn . The tempera- 
ture profile is a strong function of the crust's thermal con- 
ductivity and core's neutrino emissivity, both of which are 
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poorly constrained. We parametrize these uncertainties by 
implementing two conductivity and three core neutrino emis- 
sivity prescriptions that likely bracket their true values in 
accreting neutron stars. The thermal conductivity is a de- 
creasing function of the impurity parameter gimp = (Z 2 ) - 
(Z) 2 dltoh & Kohvamal I1993L see also Daligault & Gupta 
2009). ISchatz et all dl999lT found Q imp ~ 100 from steady- 
state nucleosynthesis calculations, altho ugh subsequent cal- 
culations suggest Qimp should be s malle r (ISchatz et al| 2003a; 
Wooslev et al. 2004a; Koi ke et alj2004tlHorowitz et alj2007l 
20091) . In additi on, fits to the qu i escent light curves 
of KS 1731-260 (IShternin et all 120071: iBrown & Cummind 
120091) and MXB 1659-29 dBrown & Cummindl2009t) require 
that Qimp ~ 1. Since both observations and molecular dy- 
namics simulations imply the crust forms an ordered lattice, 
we adopt Q lmp = 3 and 100 as the two bracketing values. 
The core neutrino emissivity, and thereby the core cooling 
rate, depends on the u nknown ultradense matter equation of 
state (for reviews, see lYakovlev & Pethicki r2004: Pag e et ail 
2006). We consider one "fast" cooling model for which 
the pion condensate process dominates and two "slow" cool- 
ing models for which either the modified Urea or nucleon- 
nucle on bremsstrahlung process dominates (see, e.g., Table 
1 of iPage et ai1l2006|): these roughly corre spond to cases 
"A," "B," and "D" of lCumming et al] ([2006, see their Table 
2). The respective core temperatures for these models are ap- 
proximately 3 x 10 7 K, 3 x 10 8 K, and 6 x 10 8 K. 
Figure [3] shows the superburst ignition column depth E; gn 

as a function of E/Endd for various neutron star models 5 . A 
1.5 MeV resonance in the 12 C + 12 C system lowers E; gn by a 
factor ks 2 and « 4 for the fiducial and maximum (w7)r val- 
ues, respectively; the lowered Ei gn values are in accord with 
the observationally inferred values for a range of realistic neu- 
tron star model parameters. Therefore, we conclude that (1) 
a strong resonance may exist at an energy w 1.5 MeV above 
the 12 C + I2 C ground state, and (2) if such a resonance ex- 
ists, it will mitigate the discrepancy between observationally 
inferred superburst ignition depths and those calculated from 
theoretical models. 

For the low-mass X-ray transient 4U 1608-522, which ex- 
hibited a superburst, the thermal quiescent luminosity con- 
strains the core temperature to be w 2.5 x 10 8 K. Fits to the su- 
perburst light curve find an ignition column Ej gn = (1 .5^. 1) x 
10 12 gem" 2 . A resonance at 1.5 MeV could make the ignition 
temperature over this range as low as (4. 1-4.8) x 10 8 K, which 
is marginally consistent with t he calculated cr ust temperature 
at the time of the superburst (Keek et al. 2007]). 

For the transient KS 1731-260, the timescale for the effec- 
tive temperature to decrease implies Qimp S 1> an d tne l° w ~ 
est observed effecti ve temperature implies that the core tem- 
peratu re is < 10 8 K ( IShternin et al.l2007l;lBrown & Cummind 
2009). Under these conditions, the temperature at £ w 
10 12 gem" 2 is unlikely to be > 3 x 10 8 K and therefore 
too cold to match the inferred ignition depth, even if the 

5 The critical E below which 12 C burns stably calculated in our global 
stability analysis is lo wer than that calculated in the one-zone model of 
ICummin g et al. 12006, co mpare to t heir Fig. 15). The reason is simple: Fol- 
lowing Cumming & Bildsten 12001), they demand that the characteristic life- 
time of a 12 C ionXc/r nuc > ? aC c- However, r nuc depends exponentially on the 
density due to plasma screening (ED, so 12 C burning occurs in a narrow 
column depth range, much like the electron captures discussed in SI2.3I Thus, 
the proper criterion is Xc/r n uc > f acc AE/E, where AE/E is similar to the 
expression given in equation ilO\ . This proper criterion gives a lower critical 
E, in accord with our results. 
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FIG. 3. — Superburst ignition column depth Ei gn as a function of the 
Eddington-scaled accretion rate S/Sejj for various model parameters. Solid 
(dashed) lines show results for models with impurity parameter Qimp = 3 
(100). "Pion," "Modified Urea," and "Bremsstrahlung" refer to the core's 
dominant neutrino emission mechanism. For a given Qimp and neutrino emis- 
sion mechanism, the three lines show results for a I2 C + l2 C reaction rate with 
no resonances (the standard rate), a hypothetical 1.5 MeV resonance with the 
fiducial strength (oj7)r = 3.4 X 10~ 8 eV, and a hypothetical 1.5 MeV reso- 
nance with an approximate maximum strength (oj7)r = 3.4 x 10~ 7 eV, from 
top to bottom. The boxes show the inferred Ej„„ and E ranges for the ma- 
jority of observed superbursts. A 1.5 MeV resonance lowers Ei™ by a factor 
m 2 and 4 for the fiducial and maximum (ut)r values, respectively. 

proposed resonance exists. Recent theoretical calculations 
of nuclear reactions in the neutron star crust suggest that 
the pyc nonuclear fusion of n eutron-rich, low-Z ions such 
as O dHorowitz et al.l |2008[) or reactions t riggered by f3- 
delayed neutron emissions dCrupta et al]|2008l) may provide 
a strong source of heating for the neutron star outer crust. In- 
deed, fits to the quiescent light curves of KS 1731-260 and 
MXB 1659-29 suggest that the heating in the outer crust 
is larger than can be acc ounted for from electron captures 
jBrown & Cumm ing 2009). Although a survey of neutron 
star models with this additional heating is outside the scope 
of this paper, we note that a strong resonance does alleviate 
the discrepancy in E; gn even if it does not entirely resolve it. 

6. SUMMARY AND DISCUSSION 

In this work, we reexamined the superburst trigger mecha- 
nism to address the discrepancy between observationally in- 
ferred superburst ignition column depths E; gn and those cal- 
cu lated in theoretical mo dels. Motivated by the suggestion 
of lKuulkers et al . (2002b) and the similarity between inferred 
ignition column depths from different sources, we first ex- 
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plored the viability of thermally unstable electron captures as 
the trigger mechanism in $2] We found that electron captures 
are always thermally stable in accreting neutron star oceans; 
thus electron captures do not trigger superbursts. We then in- 
vestigated the viability of nuclear fusion reactions other than 
12 C + 12 C. Accretion-induced nuclear reactions deplete ions 
with Z < 6 at column depths £ <C £i gn , whereas ions with 
Z > 6 fuse at S ^> Y,^. We therefore confirmed the pr oposal 
dCumming & Bildstenl200lUStrohmaver & Brownl2002l) that 
12 C + 12 C triggers superbursts. 

We then examined the I2 C + 12 C fusion rate in $4] noting 
that superburst model results would be in accord with obser- 
vations if the true fusion rate were greater than the standard 
rate by a factor > 10 2 . Two factors determine the fusion rate: 
plasma screening effects and the nuclear cross-section u(E). 
Uncertainties in, and corrections to, the plasma screening en- 
hancement to the reaction rate alter the usual enhancement by 
a factor < 2 and thus cannot resolve the discrepancy between 
superburst observations and theoretical models. However, un- 
certainties in u(E) are much larger; indeed, cr(E) is experi- 
mentally unknown at astrophysically relevant energies. We 
find that a strong resonance in the 12 C + 12 C system at an en- 
ergy near 1 .5 MeV could increase the fusion rate by a few or- 
ders of magnitude at the temperatures relevant to superbursts. 
Both theoretical optical potential models and extrapolations 
of existing experimental data suggest a resonance exists at an 
energy near 1 .5 MeV. If this is true and the resonance strength 
(w7)r is sufficiently large, it could eliminate the discrepancy 
between observationally inferred superburst ignition column 
depths and theoretical model results (see Figure |3). 

In 5JT] we outlined three fundamental problems that exist 
with superburst ignition. We address these problems below 
in the context of our results. 

(1) The results of all previous superburst models imply that 
ocean temperatures are too low for 12 C ignition at the inferred 
Sign « 10 12 gem -2 . A strong resonance near 1.5 MeV in the 
12 C + 12 C system would decrease the temperature required 
for ignition at Si gn « 10 12 gem" 2 from «6x 10 8 K to rj 5 x 
10 8 K. 

(2) Heavy-ion fusion hindrance would imply that the stan- 
dard 5-factor overestimates the true 12 C + 12 C fusion rate. 
This existence of such hind rance is currently speculative for 
both 12 C + I2 C in particular dJiang et al.l2007l) and exothermic 
fusion reactions in general ( Jiang et al. 2008; Stefan rni et al.l 
2008). Further more, the effect heavy-ion fusion hindrance 
has on resonant reactions is unknown. Therefore, it is un- 
clear whether heavy-ion fusion hindrance poses a problem for 
superburst ignition. 

(3) The 12 C yield from nucleosynthesis models is often 
lower than that required for a thermal instability. A strong res- 
onance would reduce the minimum 12 C abundance required 
for a superburst, but by only a small amount. Thus, this prob- 
lem would be attenuated but not resolved. 

Our result has implications for 12 C+ 12 C reactions in other 
contexts, namely Type la supernovae and massive stellar evo- 
lution. In addition, we have also presented a general prescrip- 
tion for understanding the screening enhancement factor in a 
multicomponent plasma. We briefly describe each of these 
topics before concluding with an outlook on future measure- 
ments. 

6.1. Implications for Type la Supernovae and Massive Stellar 

Evolution 



The fusion of 12 C is an important stage in the post-main- 
sequence evolution of a massive star, and it is the reaction 
that ignites a white dwarf and triggers a thermonuclear (Type 
la) supernova. In both systems, the competition between 
heating from the 12 C + 12 C reaction and cooling from neu- 
trino emissions determines ignition. To explore the impli- 
cations of a resonance in the reaction cross-section on these 
phenomena, we construct ignition curves (Fig.|U), defined as 
£nuc(p, T) = e„(p, T), for a 12 C- 16 plasma with X c = 0.5. We 
compute the neutrino emissivity for the pair, photo, plasma, 
and bremsstrahlung p rocesses using analytical fitting formu- 
lae (lltoh etalJ fT996). Fo r e nuc , we use the effective reac- 
tion Q-value of 9.0 MeV (IChamulak et alJ l2008). which in- 
cludes heating from both the p- and a- branches and sub- 
sequent reactions; the ignition curve is insensitive to the 
choice of Q. Three curves are plotted in Fig. H] for differ- 
ent c hoices of the 12 C+ 12 C ra te: the standard nonresonant 
rate dCaughlan & Fowled 1 19881 dotted line), a resonance at 
E R = 1 .5 MeV with our fiducial strength (tt>7) R = 3 .4 x 1 0" 8 eV 
(solid line), and a resonance at Er = 1 .5 MeV with our maxi- 
mum strength (ujj)r = 3.4 x 10~ 7 eV (dashed line). 

As is evident from Figure [4] the effect of a resonance at 
£r = 1.5 MeV is minimal for the ignition of Type la su- 
pernovae, which are tho ught to ignite at cen tral densities 
> 10 9 gem" 3 (see, e. g.. IWooslev et alj|2004bl) . It is inter- 
esting to speculate that a resonance at a lower energy might 
shift the ignition curve to lower densities. This would re- 
duce the in situ neutronization of the nickel-peak material 
synthesized in the explosive burning and th e neutronization 
during the pre-explosiye con nective burning (iPiro & Bi ldsten 
2008; Chamula k et al.l 120081) : moreover, numerical simula- 
tions dRopke et al. 2006) find that a lower central density re- 
duces the growth of the turbulent flame velocity (because the 
lower gravitational acceleration decreases the growth rate of 
the Rayleigh-Taylor instability), which leads to a less vigor- 
ous explosion and a decreased production of iron-peak ele- 
ments. Although lower densities are necessary to avoid over- 
production of neutron-rich isotopes such as 54 Fe and 58 Ni 
dWooslevI \199% llwamoto etail fl999). this may be amelio- 
rated by improved electron capture rates onto pf-she ll nuclei 
(iMartfnez-Pinedo et alll200Ct iBrachwitz et alJl2000l) . More- 
over, there is observational evidence that the majority of su- 
pernovae do und ergo electron-captur es in the innermost ss 
0.2 M Q of ejecta (Mazz aliet alJl2007l) . Given the uncertain- 
ties in modeling the progenitor evolution, flame ignition, and 
explosion, and in the dependence of the i gnition density on the 
accretion history of the white dwarf (see Lesaffre et al. 2006, 
for a recent discussion), we do not think it possible to con- 
strain the existence of such a resonance from observations at 
this time, but future modeling efforts should clearly allow for 
this possibility. 

Intriguingly, the largest effect is at p < 10 5 gem" 3 , which is 
the region encountered b y post-main-sequence massive stars 
(see lWoosley et al. 2002, and references therein). Stellar evo- 
lutionary calculations, which include the shock-induced ex- 
plosive nucleosynthesis, find that a decrease in the 12 C + 12 C 
rate leads to enhan cements in 26 Al and 60 Fe abundances 
(Gasques et al. 2007). Further calculations are needed to de- 
termine whether an enhanced 12 C + 12 C rate would produce 
interesting changes in nucleosynthesis. 

6.2. Ignition in a Multicomponent Plasma 

In jj4.1.1l we showed that heavy spectator ions increase the 
plasma screening enhancement to the thermonuclear reaction 
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FIG. 4. — Locus in the temperature-density plane where e nl | C = e„, which 
defines the ignition of 12 C for stellar burning in massive stars and for ther- 
monuclear (Type la) supernovae. The composition is 12 C-' 6 with Xq = 0.5. 
The three curves, fro m top to bottom, show e nue computed with the standard 
rate iCaughlan & Fowler 1988, dotted line), with a resonance at our fiducial 
strength, (o;7)r = 3.4 X 10 s eV (solid line), and with a resonance at our max- 
imum strength, (wy)^ = 3.4 X 10~ 7 eV (dashed line). 

rate via linear mixing rule deviations. To our knowledge, this 
work is the first to show that spectator ions affect the plasma 
screening enhancement in the thermonuclear regime. Prior 
work on the effects such deviations have on the plasma screen- 
ing enhancement focused on binary ionic mixtures consisting 
only of reactants and products (e.g. 12 C and 24 Mg) but no 
spectator ions (e.g.lOgata et al .11 9931: iDe Witt et al.ll996l) . and 
DeWitt & Slatterv (2003) concluded that linear mixing rule 
deviations always decrease the plasma screening enhance- 
ment in binary ionic mixtures (eq. l36l ). However, determin- 
ing the effect of spectator ions requires analyzing a mixture 
of three or more ions, but previous applications of linear mix- 
ing ru le deviations in ternary ionic mixtures (e. g., Oga ta et all 
1993) focuse d only on phase diagrams of cry stallizing white 
dwarfs (e. g. Jlsern et al.ll 199 it ISegretainl 19961) . not on fusion 
reactions. 

In our analysis, we tacitly assumed the plasma is uniformly 
mixed. However, recent molecular dynamics simulations 
of multicomponent plasmas exhibit cluster ing of low-Z ions 
(IWiinsch et al.l 12008b iHorowitz et al.ll2009l) . which may en- 
hance the reaction rate. This is worthy of further investiga- 
tion. 

For completeness, we note that linear mixing rule devia- 
tions are much larger for Coulomb solids (DeWi tt & Slatteryl 



2003). This has two consequences for accreting neutron stars: 
(1) Screening enhancements for multicomponent plasmas in 
the pycnonuclear regime may be orders of magnitude greater 
than currently thought. This could lower the pressures at 
which pycnonuclear reactions occur in the crust and thereby 
heat the ocean to a larger extent. (2) A multicomponent 
plasma's freezing temperature is much lower than that of a 
one-componen t plasm a. This possibly explains the results of 
IHorowitz et al.l(l2007l) . whose molecular dynamics simulation 
of a multicomponent plasma froze at T » 247 rather than the 
typical T « 175 of a one-component plasma. 

6.3. Outlook for Future Measurements 

Our conclusions are contingent on the existence of a strong 
res onan ce near 1.5 MeV in the 12 C + 12 C system. As noted 
in ^4.21 resonances exist throughout the experimentally stud- 
ied energy range and are spaced at intervals of w 0.3 MeV. 
Therefore, a resonance almost certainly exists sufficiently 
near 1.5 MeV (i.e., within the Gamow window). We cannot 
predict with confidence, however, that the resonance strength 
(w7)r is sufficiently large. Indeed, the measured resonances 
at higher energies typically increase the thermally averaged 
reaction rate (av) by a factor < 10 over the nonresonant con- 
tribution, so the resonance needs to be unusually strong. Re- 
solving this issue requires experimental measurements of the 
12 C + 12 C cross-section near the Gamow peak, which requires 
lab energies of 3 MeV with the ability to measure cross- 
sections at the 0.1 pb level. Measurements at higher ener- 
gies are probably required to map out the resonance structure 
between the Gamow peak and the other available measure- 
ments of the fusion cross-section. Such measurements are 
possible in the near term with existing laboratories and are 
certainly within reach of planned underground facilities such 
as DUSEL (t hough they will require the larger DUSEL accel- 
erator option: iGorres & Wiesch er 2006). 
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APPENDIX 
DERIVATION OF AH(Q) 

In this Appendix, we derive AH(Q), the correction to the plasma screening enhancement factor due to linear mixing rule 
deviations. Consider a strongly-coupled Coulomb plasma consisting of N = J^i^i ions, where Ni is the number of ions with 
charge Z,-, x,- = Ni/N is the n umber fraction of s pecies i, and species 1 and 2 are the reactant and product of the reaction, 
respectively, so that Z2 = 2Zi . iDeWitt et al.l d 1973b found that, neglecting quantum contributions, H(0) equals the difference in 
the Coulomb free energy before and after the reaction; since a fusion reaction destroys two reactant ions and creates one product 
ion, 

H(Q) = F^ Mm -f ex f,n ' A = F ex (Ni,N2, ...,N„)-F ex (Ni-2,N2+l,---,N„). ( A 1 ) 

Writing equation (lAlb in a more general form, 

„, m F^(N l -2AN2,N2 + AN 2: ...,N n )-F^(NuN2,...,N„) 

H(0)= AN 2 ' (A2) 
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where AN2 is the number of products created. In the limit AA^2 <C N\,N2, equation ( IA2I ) simplifies to 



H(0) = 

Using the general expression for F ex (eq. 1341 ). we find 

H(0) 



8 



8 



8N 2 dNi 



= 2/ ex (r 1 )-/ ex (r 2 )+{A/ ex -PAr x } . 



(A3) 



(A4) 



where the term 2/ ex (ri)-/ ex (T2) is the well-known result for a plasma obeying the linear mixing rule (e.g.. I Janco vicil[l977h . the 
bracketed term is AH(0)/IcbT, and we have defined the operator 



V=N 



( A_ 2 A 

\8N 2 dNi 



_8_ 

8x2 dxi 



i=i 



(A5) 



Pote khin et al.l (120091) derived an accurate, analytic fitting formula for A/ ex (see their eq. [16]). Using their formula for an un- 
polarized electron background, which is appropriate for a strongly-coupled Coulomb liquid, A/ ex = A/ ex (r, (Z), (Z 2 ), (Z 5 / 2 )), 

where T = y\- *,■]?,■ a nd (Z k ) = Y^i x iZf- From equations ( IA4l i and ( IA51 > of this work and equations (12), (14), and (16) of 
iPotekhin et aLl (f2009). we find 



AH(Q) 



= Af 



din A/ 6 * ^ <91nA/ ex . . <91nA/ ex , 2x <91nA/ ex . 5/2 / 
V\nT+ - ,„, Pln(Z) + J X>ln(Z 2 ) + J V\n(Z 5 ' 2 ) 



8\nT 



where 



<91nA/ ex _ abcT b 
8\nT 

<91nA/ e 



l+aT b ' 

I /"OH 



ain(z) 

VlnT = l + (2 5/3 -2)- 



<91n(Z 2 ) 

Vln(Z k ) = l + (2*-2) 



91n(Z 5 /2) ' 

(z*>' 



(A6) 



ain(z) 

91nA/ ex 
<91n(Z 2 } 

8 In Af x 



acT" 



1.1 + 345 3 r ( 1 

(l-S) + 0Ab lnr + 



DH 



acT' 



b r 



2C DH -C LM i+ar & 

>LM acT b 



2.2S+17S 4 
3.3+1025 3 



C L 



<91n(Z 5 /2) (^DH_^LM 1+aT b 



2.2S+17S 4 

2.2 + 68<5 3 
2.25+ ns 4 



(l-S)+0.2b lnr+ 



1-b 
1 



l-b 



+ -ln(l+ fl r i ), 



■ln(l+ar 6 ), 



{1-6) 



(A7) 
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